Syntenic Cell Wall QTLs as Versatile Breeding Tools: Intraspecific Allelic Variability and Predictability of Biomass Quality Loci in Target Plant Species.

Syntenic cell wall QTLs (SQTLs) can identify genetic determinants of biomass traits in understudied species based on results from model crops. However, their effective use in plant breeding requires SQTLs to display intraspecific allelic variability and to predict causative loci in other populations/species than the ones used for SQTLs identification. In this study, genome assemblies from different accessions of Arabidopsis, rapeseed, tomato, rice, Brachypodium and maize were used to evaluate the intraspecific variability of SQTLs. In parallel, a genome-wide association study (GWAS) on cell wall quality traits was performed in miscanthus to verify the colocalization between GWAS loci and miscanthus SQTLs. Finally, an analogous approach was applied on a set of switchgrass cell wall QTLs retrieved from the literature. These analyses revealed large SQTLs intraspecific genetic variability, ranging from presence-absence gene variation to SNPs/INDELs and changes in coded proteins. Cell wall genes displaying gene dosage regulation, such as PAL and CAD, displayed presence-absence variation in Brachypodium and rapeseed, while protein INDELs were detected for the Brachypodium homologs of the rice brittle culm-like 8 locus, which may likely impact cell wall quality. Furthermore, SQTLs significantly colocalized with the miscanthus and switchgrass QTLs, with relevant cell wall genes being retained in colocalizing regions. Overall, SQTLs are useful tools to screen germplasm for relevant genes and alleles to improve biomass quality and can increase the efficiency of plant breeding in understudied biomass crops.

Genetic dissection of morphological variation in rosette leaves and leafy heads in cabbage (Brassica oleracea var. capitata).

KEY MESSAGE: Correlations between morphological traits of cabbage rosette leaves and heads were found. Genome-wide association studies of these traits identified 50 robust quantitative trait loci in multiple years. Half of these loci affect both organs. Cabbage (Brassica oleracea var. capitata) is an economically important vegetable crop cultivated worldwide. Cabbage plants go through four vegetative stages: seedling, rosette, folding and heading. Rosette leaves are the largest leaves of cabbage plants and provide most of the energy needed to produce the leafy head. To understand the relationship and the genetic basis of leaf development and leafy head formation, 308 cabbage accessions were scored for rosette leaf and head traits in three-year field trials. Significant correlations were found between morphological traits of rosette leaves and heads, namely leaf area with the head area, height and width, and leaf width with the head area and head height, when heads were harvested at a fixed number of days after sowing. Fifty robust quantitative trait loci (QTLs) for rosette leaf and head traits distributed over all nine chromosomes were identified with genome-wide association studies. All these 50 loci were identified in multiple years and generally affect multiple traits. Twenty-five of the QTL were associated with both rosette leaf and leafy head traits. We discuss thirteen candidate genes identified in these QTL that are expressed in heading leaves, with an annotation related to auxin and other phytohormones, leaf development, and leaf polarity that likely play a role in leafy head development or rosette leaf expansion.

Diversity and Agronomic Performance of Lupinus mutabilis Germplasm in European and Andean Environments.

The introduction of Lupinus mutabilis (Andean lupin) in Europe will provide a new source of protein and oil for plant-based diets and biomass for bio-based products, while contributing to the improvement of marginal soils. This study evaluates for the first time the phenotypic variability of a large panel of L. mutabilis accessions both in their native environment and over two cropping conditions in Europe (winter crop in the Mediterranean region and summer crop in North-Central Europe), paving the way for the selection of accessions adapted to specific environments. The panel of 225 accessions included both germplasm pools from the Andean region and breeding lines from Europe. Notably, we reported higher grain yield in Mediterranean winter-cropping conditions (18 g/plant) than in the native region (9 g/plant). Instead, North European summer-cropping conditions appear more suitable for biomass production (up to 2 kg/plant). The phenotypic evaluation of 16 agronomical traits revealed significant variation in the panel. Principal component analyses pointed out flowering time, yield, and architecture-related traits as the main factors explaining variation between accessions. The Peruvian material stands out among the top-yielding accessions in Europe, characterized by early lines with high grain yield (e.g., LIB065, LIB072, and LIB155). Bolivian and Ecuadorian materials appear more valuable for the selection of genotypes for Andean conditions and for biomass production in Europe. We also observed that flowering time in the different environments is influenced by temperature accumulation. Within the panel, it is possible to identify both early and late genotypes, characterized by different thermal thresholds (600°C-700°C and 1,000-1,200°C GDD, respectively). Indications on top-yielding and early/late accessions, heritability of morpho-physiological traits, and their associations with grain yield are reported and remain largely environmental specific, underlining the importance of selecting useful genetic resources for specific environments. Altogether, these results suggest that the studied panel holds the genetic potential for the adaptation of L. mutabilis to Europe and provide the basis for initiating a breeding program based on exploiting the variation described herein.

Genotypic and tissue-specific variation of Populus nigra transcriptome profiles in response to drought.

Climate change is one of the most important challenges for mankind in the far and near future. In this regard, sustainable production of woody crops on marginal land with low water availability is a major challenge to tackle. This dataset is part of an experiment, in which we exposed three genetically differentiated genotypes of Populus nigra originating from contrasting natural habitats to gradually increasing moderate drought. RNA sequencing was performed on fine roots, developing xylem and leaves of those three genotypes under control and moderate drought conditions in order to get a comprehensive dataset on the transcriptional changes at the whole plant level under water limiting conditions. This dataset has already provided insight in the transcriptional control of saccharification potential of the three Populus genotypes under drought conditions and we suggest that our data will be valuable for further in-depth analysis regarding candidate gene identification or, on a bigger scale, for meta-transcriptome analysis.

A genome wide association study to dissect the genetic architecture of agronomic traits in Andean lupin (Lupinus mutabilis).

Establishing Lupinus mutabilis as a protein and oil crop requires improved varieties adapted to EU climates. The genetic regulation of strategic breeding traits, including plant architecture, growing cycle length and yield, is unknown. This study aimed to identify associations between 16 669 single nucleotide polymorphisms (SNPs) and 9 agronomic traits on a panel of 223 L. mutabilis accessions, grown in four environments, by applying a genome wide association study (GWAS). Seven environment-specific QTLs linked to vegetative yield, plant height, pods number and flowering time, were identified as major effect QTLs, being able to capture 6 to 20% of the phenotypic variation observed in these traits. Furthermore, two QTLs across environments were identified for flowering time on chromosome 8. The genes FAF, GAMYB and LNK, regulating major pathways involved in flowering and growth habit, as well as GA30X1, BIM1, Dr1, HDA15, HAT3, interacting with these pathways in response to hormonal and environmental cues, were prosed as candidate genes. These results are pivotal to accelerate the development of L. mutabilis varieties adapted to European cropping conditions by using marker-assisted selection (MAS), as well as to provide a framework for further functional studies on plant development and phenology in this species.

Determinants of Common Mental Disorders (CMD) among adolescent girls aged 15-19 years in Indonesia: Analysis of the 2018 National Basic Health Survey Data.

Common Mental Disorders (CMD) are distress conditions which manifest themselves with anxiety, somatic, and depressive symptoms. CMD are highly prevalent in Indonesia especially among adolescents. Adolescent girls have a higher risk to develop CMD than boys. This may be related to anaemia, potentially aggravated by early onset of menstruation. This study aimed to determine the association between haemoglobin concentration and other determinants of CMD among adolescent girls in Indonesia. Data of 1,052 adolescent girls aged 15-19 years old from the Indonesian Basic Health Survey 2018 were analysed. CMD was measured using the Self Reporting Questionnaire (SRQ-20). Principal Component Analysis of main determinants was applied and resulting principal components were investigated as risk factors for CMD. The prevalence of CMD among the study population was 16.5%. Anaemia and Age at Menarche (AAM) were not associated with CMD. Three principal components were significantly associated with higher CMD score: 1) higher parental education, better employment of the father, and living in an urban area (ß: 0.16, 95%-CI: 0.02; 0.30); 2) higher consumption of salty foods, high-fat foods, and soft drinks (ß: 0.23, 95%-CI: 0.05; 0.40); and 3) having asthma, smoking, and a higher haemoglobin concentration (ß: 1.74, 95%-CI: 1.59; 1.89). The strongest clustered associates of CMD among adolescent girls in Indonesia were asthma, smoking status, and higher haemoglobin concentration, whereas anaemia and AAM were not associated. Causality of smoking and diet to CMD could not be disentangled in this cross-sectional study. Our findings imply that adolescent girls who have asthma and smoke, as well as those having parents with higher education and secured occupation, are more likely to have mental disorders.

Common bean SNP alleles and candidate genes affecting photosynthesis under contrasting water regimes.

Water deficit is a major worldwide constraint to common bean (Phaseolus vulgaris L.) production, being photosynthesis one of the most affected physiological processes. To gain insights into the genetic basis of the photosynthetic response of common bean under water-limited conditions, a collection of 158 Portuguese accessions was grown under both well-watered and water-deficit regimes. Leaf gas-exchange parameters were measured and photosynthetic pigments quantified. The same collection was genotyped using SNP arrays, and SNP-trait associations tested considering a linear mixed model accounting for the genetic relatedness among accessions. A total of 133 SNP-trait associations were identified for net CO2 assimilation rate, transpiration rate, stomatal conductance, and chlorophylls a and b, carotenes, and xanthophyll contents. Ninety of these associations were detected under water-deficit and 43 under well-watered conditions, with only two associations common to both treatments. Identified candidate genes revealed that stomatal regulation, protein translocation across membranes, redox mechanisms, hormone, and osmotic stress signaling were the most relevant processes involved in common bean response to water-limited conditions. These candidates are now preferential targets for common bean water-deficit-tolerance breeding. Additionally, new sources of water-deficit tolerance of Andean, Mesoamerican, and admixed origin were detected as accessions valuable for breeding, and not yet explored.

Genetic Architecture of Flowering Time and Sex Determination in Hemp (Cannabis sativa L.): A Genome-Wide Association Study.

Flowering time and sex determination in hemp (Cannabis sativa L.) strongly influence fiber quality and seed production of this crop. The control of these traits is paramount for the breeding of new cultivars. Yet, little is known about the genetics underlying such complex traits and a better understanding requires in depth knowledge of the molecular mechanisms responsible for these traits. In this report, the genetic architecture of flowering time and sex determination in hemp was studied using a Genome-Wide Association Studies (GWAS) approach. Association studies were performed on a panel of 123 hemp accessions, tested in three contrasting environments, using a set of 600 K SNP markers. Altogether, eight QTLs were identified across environments; six for flowering time traits and two for sex determination. These QTLs covered genomic regions with 33 transcripts predicted to be involved in flowering and sex determination as well as a microRNA, miR156. Genes related to perception and transduction of light and transcription factors well-known to regulate flowering were identified in QTLs for flowering time traits. Transcription factors and genes involved in regulating the balance of phytohormones, specially auxins and gibberellic acid, were identified in QTLs for sex determination. Sex determination QTLs were associated with the development of male flowers in female plants and thus with the stability of sex determination in monecious plants. The present study elucidates relevant knowledge on the genetic mechanisms of flowering and sex determination traits in hemp, and provides new tools for hemp breeding.

Elucidating the Genetic Architecture of Fiber Quality in Hemp (Cannabis sativa L.) Using a Genome-Wide Association Study.

Hemp (Cannabis sativa L.) is a bast-fiber crop with a great potential in the emerging bio-based economy. Yet, hemp breeding for fiber quality is restricted and that is mainly due to the limited knowledge of the genetic architecture of its fiber quality. A panel of 123 hemp accessions, with large phenotypic variability, was used to study the genetic basis of seven cell wall and bast fiber traits relevant to fiber quality. These traits showed large genetic variance components and high values of broad sense heritability in this hemp panel, as concluded from the phenotypic evaluation across three test locations with contrasting environments. The hemp panel was genotyped using restriction site associated DNA sequencing (RAD-seq). Subsequently, a large set (> 600,000) of selected genome-wide single nucleotide polymorphism (SNP) markers was used for a genome-wide association study (GWAS) approach to get insights into quantitative trait loci (QTLs) controlling fiber quality traits. In absence of a complete hemp genome sequence, identification of QTLs was based on the following characteristics: (i) association level to traits, (ii) fraction of explained trait variance, (iii) collinearity between QTLs, and (iv) detection across different environments. Using this approach, 16 QTLs were identified across locations for different fiber quality traits, including contents of glucose, glucuronic acid, mannose, xylose, lignin, and bast fiber content. Among them, six were found across the three environments. The genetic markers composing the QTLs that are common across locations are valuable tools to develop novel genotypes of hemp with improved fiber quality. Underneath the QTLs, 12 candidate genes were identified which are likely to be involved in the biosynthesis and modification of monosaccharides, polysaccharides, and lignin. These candidate genes were suggested to play an important role in determining fiber quality in hemp. This study provides new insights into the genetic architecture of fiber traits, identifies QTLs and candidate genes that form the basis for molecular breeding for high fiber quality hemp cultivars.

Genetic Variability of Morphological, Flowering, and Biomass Quality Traits in Hemp (Cannabis sativa L.).

Hemp (Cannabis sativa L.) is a bast-fiber crop well-known for the great potential to produce sustainable fibers. Nevertheless, hemp fiber quality is a complex trait, and little is known about the phenotypic variability and heritability of fiber quality traits in hemp. The aim of this study is to gain insights into the variability in fiber quality within the hemp germplasm and to estimate the genetic components, environmental components, and genotype-by-environment (G×E) interactions on fiber quality traits in hemp. To investigate these parameters, a panel of 123 hemp accessions was phenotyped for 28 traits relevant to fiber quality at three locations in Europe, corresponding to climates of northern, central, and southern Europe. In general, hemp cultivated in northern latitudes showed a larger plant vigor while earlier flowering was characteristic of plants cultivated in southern latitudes. Extensive variability between accessions was observed for all traits. Most cell wall components (contents of monosaccharides derived from cellulose and hemicellulose; and lignin content), bast fiber content, and flowering traits revealed large genetic components with low G×E interactions and high broad-sense heritability values, making these traits suitable to maximize the genetic gains of fiber quality. In contrast, contents of pectin-related monosaccharides, most agronomic traits, and several fiber traits (fineness and decortication efficiency) showed low genetic components with large G×E interactions affecting the rankings across locations. These results suggest that pectin, agronomic traits, and fiber traits are unsuitable targets in breeding programs of hemp, as their large G×E interactions might lead to unexpected phenotypes in untested locations. Furthermore, all environmental effects on the 28 traits were statistically significant, suggesting a strong adaptive behavior of fiber quality in hemp to specific environments. The high variability in fiber quality observed in the hemp panel, the broad range in heritability, and adaptability among all traits prescribe positive prospects for the development of new hemp cultivars of excellent fiber quality.

Distribution of P1(D1) wart disease resistance in potato germplasm and GWAS identification of haplotype-specific SNP markers.

KEY MESSAGE: A Genome-Wide Association Study using 330 commercial potato varieties identified haplotype specific SNP markers associated with pathotype 1(D1) wart disease resistance. Synchytrium endobioticum is a soilborne obligate biotrophic fungus responsible for wart disease. Growing resistant varieties is the most effective way to manage the disease. This paper addresses the challenge to apply molecular markers in potato breeding. Although markers linked to Sen1 were published before, the identification of haplotype-specific single-nucleotide polymorphisms may result in marker assays with high diagnostic value. To identify hs-SNP markers, we performed a genome-wide association study (GWAS) in a panel of 330 potato varieties representative of the commercial potato gene pool. SNP markers significantly associated with pathotype 1 resistance were identified on chromosome 11, at the position of the previously identified Sen1 locus. Haplotype specificity of the SNP markers was examined through the analysis of false positives and false negatives and validated in two independent full-sib populations. This paper illustrates why it is not always feasible to design markers without false positives and false negatives for marker-assisted selection. In the case of Sen1, founders could not be traced because of a lack of identity by descent and because of the decay of linkage disequilibrium between Sen1 and flanking SNP markers. Sen1 appeared to be the main source of pathotype 1 resistance in potato varieties, but it does not explain all the resistance observed. Recombination and introgression breeding may have introduced new, albeit rare haplotypes involved in pathotype 1 resistance. The GWAS approach, in such case, is instrumental to identify SNPs with the best possible diagnostic value for marker-assisted breeding.

Assessment of heterosis in two Arabidopsis thaliana common-reference mapping populations.

Hybrid vigour, or heterosis, has been of tremendous importance in agriculture for the improvement of both crops and livestock. Notwithstanding large efforts to study the phenomenon of heterosis in the last decades, the identification of common molecular mechanisms underlying hybrid vigour remain rare. Here, we conducted a systematic survey of the degree of heterosis in Arabidopsis thaliana hybrids. For this purpose, two overlapping Arabidopsis hybrid populations were generated by crossing a large collection of naturally occurring accessions to two common reference lines. In these Arabidopsis hybrid populations the range of heterosis for several developmental and yield related traits was examined, and the relationship between them was studied. The traits under study were projected leaf area at 17 days after sowing, flowering time, height of the main inflorescence, number of side branches from the main stem or from the rosette base, total seed yield, seed weight, seed size and the estimated number of seeds per plant. Predominantly positive heterosis was observed for leaf area and height of the main inflorescence, whereas mainly negative heterosis was observed for rosette branching. For the other traits both positive and negative heterosis was observed in roughly equal amounts. For flowering time and seed size only low levels of heterosis were detected. In general the observed heterosis levels were highly trait specific. Furthermore, no correlation was observed between heterosis levels and the genetic distance between the parental lines. Since all selected lines were a part of the Arabidopsis genome wide association (GWA) mapping panel, a genetic mapping approach was applied to identify possible regions harbouring genetic factors causal for heterosis, with separate calculations for additive and dominance effects. Our study showed that the genetic mechanisms underlying heterosis were highly trait specific in our hybrid populations and greatly depended on the genetic background, confirming the elusive character of heterosis.

Physical mapping of QTL for tuber yield, starch content and starch yield in tetraploid potato (Solanum tuberosum L.) by means of genome wide genotyping by sequencing and the 8.3 K SolCAP SNP array.

BACKGROUND: Tuber yield and starch content of the cultivated potato are complex traits of decisive importance for breeding improved varieties. Natural variation of tuber yield and starch content depends on the environment and on multiple, mostly unknown genetic factors. Dissection and molecular identification of the genes and their natural allelic variants controlling these complex traits will lead to the development of diagnostic DNA-based markers, by which precision and efficiency of selection can be increased (precision breeding). RESULTS: Three case-control populations were assembled from tetraploid potato cultivars based on maximizing the differences between high and low tuber yield (TY), starch content (TSC) and starch yield (TSY, arithmetic product of TY and TSC). The case-control populations were genotyped by restriction-site associated DNA sequencing (RADseq) and the 8.3 k SolCAP SNP genotyping array. The allele frequencies of single nucleotide polymorphisms (SNPs) were compared between cases and controls. RADseq identified, depending on data filtering criteria, between 6664 and 450 genes with one or more differential SNPs for one, two or all three traits. Differential SNPs in 275 genes were detected using the SolCAP array. A genome wide association study using the SolCAP array on an independent, unselected population identified SNPs associated with tuber starch content in 117 genes. Physical mapping of the genes containing differential or associated SNPs, and comparisons between the two genome wide genotyping methods and two different populations identified genome segments on all twelve potato chromosomes harboring one or more quantitative trait loci (QTL) for TY, TSC and TSY. CONCLUSIONS: Several hundred genes control tuber yield and starch content in potato. They are unequally distributed on all potato chromosomes, forming clusters between 0.5-4 Mbp width. The largest fraction of these genes had unknown function, followed by genes with putative signalling and regulatory functions. The genetic control of tuber yield and starch content is interlinked. Most differential SNPs affecting both traits had antagonistic effects: The allele increasing TY decreased TSC and vice versa. Exceptions were 89 SNP alleles which had synergistic effects on TY, TSC and TSY. These and the corresponding genes are primary targets for developing diagnostic markers.

Targeted and Untargeted Approaches Unravel Novel Candidate Genes and Diagnostic SNPs for Quantitative Resistance of the Potato (Solanum tuberosum L.) to Phytophthora infestans Causing the Late Blight Disease.

The oomycete Phytophthora infestans causes late blight of potato, which can completely destroy the crop. Therefore, for the past 160 years, late blight has been the most important potato disease worldwide. The identification of cultivars with high and durable field resistance to P. infestans is an objective of most potato breeding programs. This type of resistance is polygenic and therefore quantitative. Its evaluation requires multi-year and location trials. Furthermore, quantitative resistance to late blight correlates with late plant maturity, a negative agricultural trait. Knowledge of the molecular genetic basis of quantitative resistance to late blight not compromised by late maturity is very limited. It is however essential for developing diagnostic DNA markers that facilitate the efficient combination of superior resistance alleles in improved cultivars. We used association genetics in a population of 184 tetraploid potato cultivars in order to identify single nucleotide polymorphisms (SNPs) that are associated with maturity corrected resistance (MCR) to late blight. The population was genotyped for almost 9000 SNPs from three different sources. The first source was candidate genes specifically selected for their function in the jasmonate pathway. The second source was novel candidate genes selected based on comparative transcript profiling (RNA-Seq) of groups of genotypes with contrasting levels of quantitative resistance to P. infestans. The third source was the first generation 8.3k SolCAP SNP genotyping array available in potato for genome wide association studies (GWAS). Twenty seven SNPs from all three sources showed robust association with MCR. Some of those were located in genes that are strong candidates for directly controlling quantitative resistance, based on functional annotation. Most important were: a lipoxygenase (jasmonate pathway), a 3-hydroxy-3-methylglutaryl coenzyme A reductase (mevalonate pathway), a P450 protein (terpene biosynthesis), a transcription factor and a homolog of a major gene for resistance to P. infestans from the wild potato species Solanum venturii. The candidate gene approach and GWAS complemented each other as they identified different genes. The results of this study provide new insight in the molecular genetic basis of quantitative resistance in potato and a toolbox of diagnostic SNP markers for breeding applications.

Population structure and genome-wide association analysis for frost tolerance in oat using continuous SNP array signal intensity ratios.

KEY MESSAGE: Infinium SNP data analysed as continuous intensity ratios enabled associating genotypic and phenotypic data from heterogeneous oat samples, showing that association mapping for frost tolerance is a feasible option. Oat is sensitive to freezing temperatures, which restricts the cultivation of fall-sown or winter oats to regions with milder winters. Fall-sown oats have a longer growth cycle, mature earlier, and have a higher productivity than spring-sown oats, therefore improving frost tolerance is an important goal in oat breeding. Our aim was to test the effectiveness of a Genome-Wide Association Study (GWAS) for mapping QTLs related to frost tolerance, using an approach that tolerates continuously distributed signals from SNPs in bulked samples from heterogeneous accessions. A collection of 138 European oat accessions, including landraces, old and modern varieties from 27 countries was genotyped using the Infinium 6K SNP array. The SNP data were analyzed as continuous intensity ratios, rather than converting them into discrete values by genotype calling. PCA and Ward's clustering of genetic similarities revealed the presence of two main groups of accessions, which roughly corresponded to Continental Europe and Mediterranean/Atlantic Europe, although a total of eight subgroups can be distinguished. The accessions were phenotyped for frost tolerance under controlled conditions by measuring fluorescence quantum yield of photosystem II after a freezing stress. GWAS were performed by a linear mixed model approach, comparing different corrections for population structure. All models detected three robust QTLs, two of which co-mapped with QTLs identified earlier in bi-parental mapping populations. The approach used in the present work shows that SNP array data of heterogeneous hexaploid oat samples can be successfully used to determine genetic similarities and to map associations to quantitative phenotypic traits.

Variation in broccoli cultivar phytochemical content under organic and conventional management systems: implications in breeding for nutrition.

Organic agriculture requires cultivars that can adapt to organic crop management systems without the use of synthetic pesticides as well as genotypes with improved nutritional value. The aim of this study encompassing 16 experiments was to compare 23 broccoli cultivars for the content of phytochemicals associated with health promotion grown under organic and conventional management in spring and fall plantings in two broccoli growing regions in the US (Oregon and Maine). The phytochemicals quantified included: glucosinolates (glucoraphanin, glucobrassicin, neoglucobrassin), tocopherols (δ-, γ-, α-tocopherol) and carotenoids (lutein, zeaxanthin, ß-carotene). For glucoraphanin (17.5%) and lutein (13%), genotype was the major source of total variation; for glucobrassicin, region (36%) and the interaction of location and season (27.5%); and for neoglucobrassicin, both genotype (36.8%) and its interactions (34.4%) with season were important. For δ- and γ-tocopherols, season played the largest role in the total variation followed by location and genotype; for total carotenoids, genotype (8.41-13.03%) was the largest source of variation and its interactions with location and season. Overall, phytochemicals were not significantly influenced by management system. We observed that the cultivars with the highest concentrations of glucoraphanin had the lowest for glucobrassicin and neoglucobrassicin. The genotypes with high concentrations of glucobrassicin and neoglucobrassicin were the same cultivars and were early maturing F1 hybrids. Cultivars highest in tocopherols and carotenoids were open pollinated or early maturing F1 hybrids. We identified distinct locations and seasons where phytochemical performance was higher for each compound. Correlations among horticulture traits and phytochemicals demonstrated that glucoraphanin was negatively correlated with the carotenoids and the carotenoids were correlated with one another. Little or no association between phytochemical concentration and date of cultivar release was observed, suggesting that modern breeding has not negatively influenced the level of tested compounds. We found no significant differences among cultivars from different seed companies.

Comparative methods for association studies: a case study on metabolite variation in a Brassica rapa core collection.

BACKGROUND: Association mapping is a statistical approach combining phenotypic traits and genetic diversity in natural populations with the goal of correlating the variation present at phenotypic and allelic levels. It is essential to separate the true effect of genetic variation from other confounding factors, such as adaptation to different uses and geographical locations. The rapid availability of large datasets makes it necessary to explore statistical methods that can be computationally less intensive and more flexible for data exploration. METHODOLOGY/PRINCIPAL FINDINGS: A core collection of 168 Brassica rapa accessions of different morphotypes and origins was explored to find genetic association between markers and metabolites: tocopherols, carotenoids, chlorophylls and folate. A widely used linear model with modifications to account for population structure and kinship was followed for association mapping. In addition, a machine learning algorithm called Random Forest (RF) was used as a comparison. Comparison of results across methods resulted in the selection of a set of significant markers as promising candidates for further work. This set of markers associated to the metabolites can potentially be applied for the selection of genotypes with elevated levels of these metabolites. CONCLUSIONS/SIGNIFICANCE: The incorporation of the kinship correction into the association model did not reduce the number of significantly associated markers. However incorporation of the STRUCTURE correction (Q matrix) in the linear regression model greatly reduced the number of significantly associated markers. Additionally, our results demonstrate that RF is an interesting complementary method with added value in association studies in plants, which is illustrated by the overlap in markers identified using RF and a linear mixed model with correction for kinship and population structure. Several markers that were selected in RF and in the models with correction for kinship, but not for population structure, were also identified as QTLs in two bi-parental DH populations.

Analysis of natural allelic variation in Arabidopsis using a multiparent recombinant inbred line population.

To exploit the diversity in Arabidopsis thaliana, eight founder accessions were crossed to produce six recombinant inbred line (RIL) subpopulations, together called an Arabidopsis multiparent RIL (AMPRIL) population. Founders were crossed pairwise to produce four F1 hybrids. These F1s were crossed according to a diallel scheme. The resulting offspring was then selfed for three generations. The F4 generation was genotyped with SNP and microsatellite markers. Data for flowering time and leaf morphology traits were determined in the F5 generation. Quantitative trait locus (QTL) analysis for these traits was performed using especially developed mixed-model methodology, allowing tests for QTL main effects, QTL by background interactions, and QTL by QTL interactions. Because RILs were genotyped in the F4 generation and phenotyped in the F5 generation, residual heterozygosity could be used to confirm and fine-map a number of the QTLs in the selfed progeny of lines containing such heterozygosity. The AMPRIL population is an attractive resource for the study of complex traits.

Statistical epistasis between candidate gene alleles for complex tuber traits in an association mapping population of tetraploid potato.

Association mapping using DNA-based markers is a novel tool in plant genetics for the analysis of complex traits. Potato tuber yield, starch content, starch yield and chip color are complex traits of agronomic relevance, for which carbohydrate metabolism plays an important role. At the functional level, the genes and biochemical pathways involved in carbohydrate metabolism are among the best studied in plants. Quantitative traits such as tuber starch and sugar content are therefore models for association genetics in potato based on candidate genes. In an association mapping experiment conducted with a population of 243 tetraploid potato varieties and breeding clones, we previously identified associations between individual candidate gene alleles and tuber starch content, starch yield and chip quality. In the present paper, we tested 190 DNA markers at 36 loci scored in the same association mapping population for pairwise statistical epistatic interactions. Fifty marker pairs were associated mainly with tuber starch content and/or starch yield, at a cut-off value of q ≤ 0.20 for the experiment-wide false discovery rate (FDR). Thirteen marker pairs had an FDR of q ≤ 0.10. Alleles at loci encoding ribulose-bisphosphate carboxylase/oxygenase activase (Rca), sucrose phosphate synthase (Sps) and vacuolar invertase (Pain1) were most frequently involved in statistical epistatic interactions. The largest effect on tuber starch content and starch yield was observed for the paired alleles Pain1-8c and Rca-1a, explaining 9 and 10% of the total variance, respectively. The combination of these two alleles increased the means of tuber starch content and starch yield. Biological models to explain the observed statistical epistatic interactions are discussed.

Natural DNA variation at candidate loci is associated with potato chip color, tuber starch content, yield and starch yield.

Complex characters of plants such as starch and sugar content of seeds, fruits, tubers and roots are controlled by multiple genetic and environmental factors. Understanding their molecular basis will facilitate diagnosis and combination of superior alleles in crop improvement programs ("precision breeding"). Association genetics based on candidate genes is one approach toward this goal. Tetraploid potato varieties and breeding clones related by descent were evaluated for 2 years for chip quality before and after cold storage, tuber starch content, yield and starch yield. Chip quality is inversely correlated with tuber sugar content. A total of 36 loci on 11 potato chromosomes were evaluated for natural DNA variation in 243 individuals. These loci included microsatellites and genes coding for enzymes that function in carbohydrate metabolism or transport (candidate loci). The markers were used to analyze population structure and were tested for association with the tuber quality traits. Highly significant and robust associations of markers with 1-4 traits were identified. Most frequent were associations with chip quality and tuber starch content. Alleles increasing tuber starch content improved chip quality and vice versa. With two exceptions, the most significant and robust associations (q < 0.01) were observed with DNA variants in genes encoding enzymes that function in starch and sugar metabolism or transport. Comparing linkage and linkage disequilibrium between loci provided evidence for the existence of large haplotype blocks in the breeding materials analyzed.